Energy and centrality dependence of p and p production and the A/p ratio 
in Pb+Pb collisions between 20A GeV and 158A GeV 
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■ The transverse mass (mt) distributions for antiprotons are measured at midrapidity for minimum 

bias Pb+Pb collisions at 158,4 GeV and for central Pb+Pb collisions at 20, 30, 40 and 80 A GeV 
beam energies in the fixed target experiment NA49 at the CERN SPS. The rapidity density dn/dy, 
inverse slope parameter T and mean transverse mass {m t ) derived from the m t distributions are 
studied as a function of the incident energy and the collision centrality and compared to the relevant 
data on proton production. The shapes of the mt distributions of p and p are very similar. The 
ratios of the particle yields, p/p and A/p, are also analysed. The p/p ratio exhibits an increase with 
decreasing centrality and a steep rise with increasing beam energy. The A/p ratio increases beyond 
unity with decreasing beam energy. 
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INTRODUCTION 

Copious antibaryon as well as strangeness production 
in relativistic heavy ion collisions relative to the corre- 
sponding yields observed in elementary hadronic interac- 
tions have been suggested as signatures of the QCD phase 
transition to a deconfined partonic state, the quark-gluon 
plasma 0, Q • The enhancement was expected to arise 



from gluon fragmentation into quark-antiquark pairs 
which is believed to have a significantly lower threshold 
than baryon-antibaryon and strange-antistrange hadron 
pair production channels. More generally, such enhance- 
ments should be a consequence of the creation of a large 
volume of high energy density matter, uniquely charac- 
teristic of central relativistic nuclear collisions |2| • 

Meanwhile, in the net-baryon rich systems, significant 
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annihilation losses could occur before antibaryons escape 
the collision volume. Thus a systematic study of antipro- 
ton production has been proposed as an indirect way of 
measuring the baryon density created at the instant of 
hadron formation |g. The latter is related to the degree 
of baryon stopping achieved during the early interpene- 
tration phase of the collision. It has also been pointed 
out that antiprotons can be regenerated in multi-meson 
reactions during the final hadron gas stage Experi- 
mentally, the interplay between baryon stopping, baryon 
pair production and antibaryon annihilation can be as- 
sessed through the measurement of antiproton and pro- 
ton yields, and the yield ratios. This is the subject of the 
study presented in this paper for the CERN SPS energy 
range. 

As no antibaryons are contained in the initial nuclear 
projectiles their yields and spectra are determined pre- 
dominantly by processes occurring subsequent to the 
primordial baryon stopping mechanisms. Of particu- 
lar interest is the ratio of antilambda (A) to antiproton 
(p) which, intuitively, reflects the relative abundances 
of anti-strange to anti-light quarks at the stage of for- 
mation of the finally observed antibaryons. The A/p 
production ratio amounts to about 0.25 in elementary 
hadron collisions at SPS energy but it has been found 
to significantly exceed unity in previous studies of cen- 
tral nucleus-nucleus collisions at AGS energies 0, S • 
From top SPS energy ( 1 /s^T=l7.3 GeV) to RHIC en- 
ergy (V s ivjv = 130 GeV) this ratio is known to stay at 
about unity 0, 0] , a value also obtained as an upper 
limit (reflecting the general observation on strangeness 
enhancement in central A+A collisions) in recent the- 
oretical studies emplo yin g statistical, and microscopic 
transport models 0, 0, Q, 0, 0|. As, in particu- 
lar, the latter dynamical models attempt to incorporate 
both the effects of primordial net-baryon stopping and 
of antibaryon annihilation (subsequent to the initial an- 
tibaryon formation phase) [lallflUjlllii the observation 
at top AGS energies, of A/p significantly exceeding unity, 
remains a puzzle. This effect might be related to the re- 
cent observation of a steep maximum of the K + /tt + ratio, 
occuring in the vicinity of top AGS to lowest SPS energy 
[lil EfT], Similar to the A/p ratio, the K + /tt + ratio also 
refers to the anti-strange to anti-light quark abundance 
ratio, prevaling at hadron formation. 

The paper presents an extension of the previously 
known p/p production ratio to the lower CERN SPS ener- 
gies. New measurements of antiproton yields in centrality 
selected Pb+Pb collisions at beam energies of 158^4 GeV 
(^s Niv = 17.3 GeV) and in central Pb+Pb collisions at 
20, 30, 40 and 80 A GeV ( v /s^7=6.3, 7.6, 8.7 and 12.3 
GeV, respectively) were performed bridging the gap be- 
tween the data at top AGS and SPS energies. Previously 
published proton yield measurements of NA49 at 40, 80 
and 158 A GeV [2l| are extended to 20 and 30 A GeV 
beam energies. The A/p ratios are then obtained using 
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EXPERIMENTAL METHOD 

NA49 is a fixed target exp eriment at CERN with a 
large acceptance detector [2^ using external SPS beams 
of nuclei and hadrons. The produced charged parti- 
cles were detected in four large volume Time Projection 
Chambers (TPCs). Two of them (VTPCs) are located 
inside the magnetic field, the two others (MTPCs) down- 
stream of the magnets on either side of the beam line. 
The TPCs provide precise tracking and particle identifi- 
cation in a wide range of phase space based on the mea- 
sured momentum and the specific energy loss dE/dx in 
the TPC gas with about 4% resolution. 

Two TOF detector walls containing 891 scintillator 
pixels each are situated symmetrically behind the TPCs 
on both sides of the beam. The average overall time res- 
olution of these detectors is 60 — 70 ps. In the present 
analysis, the identification of antiprotons was primarily 
accomplished by the TOF measurement, the dE/dx in- 
formation from the large TPCs being employed to reduce 
the background of charged pions and kaons in the mass 
spectrum. To demonstrate the identification capability, 
a sample of the particle mass spectra obtained from mea- 
sured momenta and time-of-flight are shown in Fig. 1. 

On-line event characterization and triggering is accom- 
plished by beam definition detectors located in the beam 
line upstream of the target and interaction counters, and 
a calorimeter downstream of the target. The data sam- 
ples were recorded with two trigger settings, providing 
the selection of central and minimum bias events. 

Central Pb+Pb collisions were selected by request- 
ing the energy deposited in the projectile fragmentation 
region to be lower than a given threshold. This was 
achieved with a Zero Degree Calorimeter (ZDC) located 
downstream of the detector measuring the energy Ez dc 
of the remaining projectile fragments and spectator pro- 
tons and neutrons. The upper limit on the energy in the 
ZDC was set to accept the 12% most central events at 
158 A GeV and 7% at 20, 30, 40 and 80 A GeV from all 
inelastic Pb+Pb collisions. For the top SPS energy, the 
5% most central interactions were selected offline. 

A Gas Cherenkov counter provides a minimum bias 
trigger for Pb+Pb collisions. It is placed in the gas region 
immediately behind the target to veto non-interacting 
projectiles. Triggering is accomplished by placing an up- 
per threshold on the signal from this Cherenkov counter 
in coincidence with a valid signal from the beam detec- 
tors. Additionally, off-line cuts were made on the position 
of the fitted primary vertex along the beam direction to 
minimize the fraction of non-target background events. 
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ons in the geometrical overlap volume of the colliding 
nuclei. The values of N woun d as a function of impact pa- 
rameter b and inelastic cross section range are obtained 
from Glauber model calculations [25J. The resulting av- 
erages (N woun d} and (6) for the six inelastic cross section 
intervals of this analysis are shown in Table I. The value 
of <7i ne i is not known precisely. An estimated uncertainty 
of 5% leads to a systematic uncertainty in (N woun d) of 
1% ( 10% ) for central (peripheral) collisions. Note that 
bin 6 suffers from a trigger bias, i.e. the minimum bias 
trigger only accepts about 35 % of the events in this cen- 
trality range. The listed values of (N woun d) and (N part ) 
are not averages for the 43 - 100% cross section interval, 
but refer to those events accepted by the trigger. Since 
the quantity N woun d is not a direct observable, previous 
publications of NA49 also employed the number of par- 
ticipants N part , which counts all the nucleons from the 
incident nuclei involved in any hadronic reaction (includ- 
ing those hit by a scattered or produced particle). N part 
can be estimated from the baryons and antibaryons mea- 
sured by the detector or from the spectator energy de- 
posited in the ZDC [H HI . For comparison (N part ) 
is also listed in Table I; the values differ significantly from 
{N woun d} for peripheral collisions. 

TABLE I. The six centrality classes used in the anal- 
ysis for 158A GeV Pb+Pb collisions, listing the covered 
range in fraction of the total cross section <Ji ne i , the mean 
numbers of participating (N part ) and wounded (N woun d) 
nucleons, and impact parameter (b) for the corresponding 
cross sections. The standard deviations of the distribu- 
tions of N woun d and b for the respective centrality bins 
are also quoted. 



FIG. 1: Mass squared histograms for TOF-identified kaons 
and antiprotons in Pb+Pb reactions at 20, 30, 40, 80 and 158 
A GeV beam energies. dE/dx cuts have been applied to the 
data to reject most of the pions. The particle momenta range 
from 3 to 10 GeV/c. 



DATA ANALYSIS AND RESULTS 

For the analysis of 158A GeV Pb+Pb collisions the 
data samples of 450,000 central and 400,000 minimum 
bias events were used. In order to study the centrality 
dependence the data were segregated into six centrality 
classes, which were selected by subdividing the range of 
the total energy measured in the ZDC. For each central- 
ity bin the cross section is determined experimentally and 
corresponds to an interval in the inelstic Pb+Pb collision 
cross section. The resulting intervals, assuming a total 
inelastic cross section Gi ne i =7.15 barns, are listed in Ta- 
ble I. The initial geometry of nucleus- nucleus collisions is 
usually characterized by the number of wounded nucleons 
N W ound which essentially represent the number of nucle- 



Centrality 


Fraction 


(Npart) 


(-^V- wound) 


(b) 


bin 


of a ine i{%) 






(fm) 


1 


0-5 


366 


357 ±22 


2.3 ±0.9 


2 


5-12 


309 


288 ± 28 


4.5 ±0.8 


3 


12-23 


242 


211 + 30 


6.4 ±0.8 


4 


23-33 


178 


146 ±25 


8.1 ±0.8 


5 


33-43 


132 


99 ±22 


9.5 ±0.8 


G 


43-100 


85 


42 ± 16 


11.8 ± 1.0 



For the reconstruction of charged tracks, a global 
tracking scheme was used which combined all track seg- 
ments reconstructed in the TPCs that belong to the same 
particle. For each event the primary vertex was de- 
termined from the intersection of reconstructed tracks. 
Events in which no primary vertex was found were re- 
jected. Event vertex as well as track quality cuts were 
applied in order to select the events for further analysis. 
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The event vertex had to lie within ±1 cm of the target 
foil position. The global track was required to comprise 
track segments in the MTPC and at least one of the VT- 
PCs. The reconstructed tracks having momentum and 
dE/dx information were extrapolated to the TOF detec- 
tor wall and assigned a mass squared value m 2 derived 
from the momentum, flight path and time-of-flight mea- 
surements. Furthermore, cuts were applied to eliminate 
tracks which impinge close to the edges of the scintillator 
tiles or which point to tiles hit by more than one track or 
which show a signal from event-correlated 7-conversion 
pairs in the scintillator. The corresponding inefficiency 
was determined experimentally from the TPC tracking 
data (edge and double hit) and the charge measurements 
in the TOF scintillators (7 conversion) and amounts to 
25% on average, with a maximum of 30% in the region of 
the TOF wall closest to the beam. The relevant correc- 
tions have been applied to the measured particle yields. 

Next, antiprotons were selected and their transverse 
mass nit distributions reconstructed. The tracks were 
subjected to identification cuts applied to the measured 
dE/dx and m 2 values in order to simultaneously maxi- 
mize the number of antiprotons and minimize the back- 
ground stemming mostly from pions and kaons. Correc- 
tion factors for the p yield due to the cuts were estimated 
from the experimental dE/dx and m 2 distributions using 
parametrized descriptions obtained from a fit. The mass 
squared distributions were described by a sum of a Gaus- 
sian distribution for the p signal and a background with 
two components: an exponentially falling contribution 
from the tails of K~ and n~ bands, and a flat random 
distribution due to misidentified particles. The energy 
loss dE/dx is well described by a sum of Gaussian dis- 
tributions. The fits were performed for each (p t ,p)-bin 
of width Ap f =0.2 GeV/c and Ap=1.0 GeV/c in the full 
kinematic range of antiprotons accepted by the TOF de- 
tector: < p t < 1.7 GeV/c and 3.0 < p < 10.0 GeV/c. 
The correction factors for the antiproton loss and back- 
ground contamination due to the cuts were then deter- 
mined separately for each (p t ,p)-bin. They were negli- 
gible in the momentum range up to p = 6.0 GeV/c and 
reached about 15% at large momenta. 

Some fraction of the measured antiprotons are the 
daughters of weak decays of strange antibaryons, mainly 
A, including those from electromagnetic decay of E° 
which are experimentally indistinguishable from primary 
A. This so called feed-down was evaluated from a 
GEANT-based Monte Carlo simulation of A and E de- 
cays in the NA49 detector, including detector response 
simulation and reconstruction of the charged decay prod- 
ucts. This procedure takes into account the probability 
for a secondary antiproton to be reconstructed as a track 
from the primary vertex. As input for the simulation, 
A yields and phase space distributions measured by the 
SPS experiments were employed. The A yields in each 
centrality bin were derived using the A centrality depen- 



dence obtained in NA57 [23| scaled to the A yield in cen- 
tral Pb+Pb collisions measured by NA49 22]. The rela- 
tively small fraction of antiprotons from E~ -decays was 
calculated using the RQMD model |29J simulation. The 
calculated m t -spectra of feeddown antiprotons are well 
described by an exponential function. The inverse slope 
parameter of these spectra gradually changes from 265 
MeV to 215 MeV from the most central to the peripheral 
bins, respectively. The overall contribution of feed-down 
antiprotons was found to vary from 25% in peripheral to 
35% in central collisions. These values contain also the 
antiprotons from doubly strange cascades. 

The correction factor for geometrical acceptance was 
calculated using the GEANT package for tracking par- 
ticles and dedicated NA49 software to simulate the de- 
tector response. The track reconstruction efficiency in 
the TPCs for primary p and p was determined by em- 
bedding simulated particle tracks into raw data events, 
which were then passed through the same reconstruction 
procedure as the real data. It was found to be nearly 
100% in the kinematical range covered by the TOF ac- 
ceptance. 

The discussed corrections have been applied to the 
data in each rapidity- transverse mass (y, m t )-bin. The 
transverse mass nit [nit = yPt + m2 ) spectra for an- 
tiprotons were then obtained by integrating the data 
over the measured rapidity region. The resulting p spec- 
tra as well as previously published proton spectra |2l) 
in centrality selected 158A GeV Pb+Pb collisions at 
2.4 < y < 2.8 (yem — 2.9) are shown in Fig. 2 along 
with a fit function of the form: 

d2n = C 1 e-^- m ^ T + C 2 e^ m '- m ^ T ' , (1) 
m t dm t dy 

A single exponential (first term in Eq. (1)) well de- 
scribes the data at m t — m > 0.2 GeV/c 2 and determines 
the inverse slope parameter of the spectra to be discussed 
further (hereafter denoted by T). The second exponential 
term with a slope paremeter T" « 100 MeV accounts for 
the flattening of the spectra at low m t . Its contribution 
to the total yield decreases from nearly 10% in central 
events to almost zero in the peripheral bins. The double 
exponential parameterization of Eq. (1) was used for ex- 
trapolation to the unmeasured m t region when calculat- 
ing the midrapidity yield dn / dy and the mean transverse 
mass (nit)—m. The extrapolation of the antiproton spec- 
tra ranges from 8% to 15% for the two most central bins 
and for the more peripheral bins, respectively. Table II 
summarizes the results on dn/dy, T and (m t )—ni derived 
from the p and p transverse spectra for each centrality bin 
of 158A GeV Pb+Pb collisions. 

The sources of considered systematic uncertainties in- 
clude possible errors in the determination of the efficiency 
corrections, in particular those for selection of single 
track hits in the TOF pixels as well as the corrections for 
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FIG. 2: (Color online) Transverse mass distributions for an- 
tiprotons and protons |2l| in the rapidity interval 2.4 < y < 
2.8 for the six centrality classes in 158,4 GeV Pb+Pb colli- 
sions. The circle and square points denote the data obtained 
with the trigger settings selecting minimum bias and central 
collisions, respectively. The errors are statistical. The curves 
through the measured points represent the two-exponential 
fit of Eq. (1) to the data. For clarity, the spectra are scaled 
down by a factor of 2 successively from the uppermost data. 



particle identification and background subtraction. This 
is only relevant to the antiprotons at higher momenta 
and high particle multiplicity. 

The estimate of the particle yield dn / dy represents an 
extrapolation into unmeasured regions under the explicit 
assumption of a certain shape of the mt distribution. 
However, the major systematic error arises from uncer- 



tainty in the phase space distributions of the parent hy- 
perons decaying into antiprotons. The latter was derived 
from the spread of the experimental data on A produc- 
tion used for feed-down corrections including the quoted 
systematic errors. The overall systematic errors of the 
antiproton results amount to approximately 10% for the 
particle yield and 5% for the inverse slope parameter with 
a weak dependence on centrality of the collision. 

Analyses of transverse spectra at the lower SPS en- 
ergies of 20, 30, 40 and 80 A GeV were performed in 
the same manner as for 158,4 GeV. Here, the online 
trigger was set to the 7% most central Pb+Pb colli- 
sions which corresponds to the number of wounded nu- 
cleons (N woun d)—34:5±28. The data set comprises about 
300,000 events for each energy. Measurements were done 
near midrapidity and covered the rapidity intervals of 
1.5 < y < 2.2 for 20,4 GeV (y cm = 1.88), 1.6 < y < 2.3 
for 30A GeV (y cm = 2.08), 1.9 < y < 2.3 for 40,4 
GeV (y cm = 2.22) and 2.2 < y < 2.6 for 80.4 GeV 
{Vcm = 2.57). 

TABLE II. Particle yield dn/dy, inverse slope T and 
mean transverse mass (mt) — m for antiprotons and pro- 
tons (2.4 < y < 2.8) at various centralities in 158,4 GeV 
Pb+Pb collisions. The errors are statistical. Proton re- 
sults are from fell. 





Centrality 

(% Of (Tinel) 


dn/dy 


T 
(MeV) 


(m t ) - m 
(MeV/c 2 ) 


p 


0-5 


1.66 + 0.17 


291 + 15 


384+ 19 




5-12 


1.27+0.11 


299+ 15 


393+ 16 




12-23 


1.05 + 0.08 


274 ± 22 


370 ± 45 




23-33 


0.76 + 0.06 


269 ± 29 


320 ± 35 




33-43 


0.55 + 0.05 


255 ± 28 


309 ± 32 




43-100 


0.33 + 0.04 


236 ± 26 


284 ± 28 


p 


0-5 


29.6 + 0.9 


308 + 9 


413+13 




5-12 


22.2 + 0.6 


308 + 9 


415+14 




12-23 


14.5 + 0.4 


276 + 9 


362+ 12 




23-33 


9.8 + 0.3 


273 ± 10 


355+ 12 




33-43 


5.7 + 0.2 


245 ± 10 


315+13 




43-100 


2.9 + 0.1 


216+10 


259+ 12 



The geometrical acceptance at lower energies allows to 
extend the transverse mass distributions for protons up 
to mt - m = 1.5 GeV/c 2 (p t = 2.25 GeV/c). Thus, the 
extrapolation into the unmeasured mt region is negligibly 
small. The background contamination due to misidenti- 
fication of antiprotons at lower energies is less than that 
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FIG. 3: (Color online) Transverse mass distributions for an- 
tiprotons and protons at midrapidity in central Pb+Pb col- 
lisions at 20, 30, 40, 80 and 158 A GeV. The errors are sta- 
tistical. The solid lines illustrate the two-exponential fit of 
Eq. (1) to the data. For clarity, the spectra are scaled down 
by a factor of 2 successively from the uppermost data. 

at 158 A GeV. 

Corrections due to the identification cuts applied to 
dE/dx and m 2 for selection of the p samples do not 
exceed a few percent at the largest values of p and pt- 
For the feed-down correction, the NA49 data |22ll30| on 
A and A production obtained for the same data sam- 
ples were used. The p feed-down contribution at lower 
energies became larger as compared to 158A GeV and 
gradually increases with decreasing beam energy reaching 



nearly 50% at 20 and 30 A GeV. This increases the over- 
all systematic errors towards lower energies approaching 
15% for the antiproton yield at 20A GeV. The feed-down 
correction for p at low energies amounts to approximately 
15% which is rather close to that for 158A GeV. 

Fig. 3 depicts the p spectra for central Pb+Pb colli- 
sions at all five energies together with fits by the double 
exponential function Eq. (1). A deviation from a single 
exponential shape is clearly seen. The slope parameter 
T" w 100 MeV in the second term of Eq.(l) is almost the 
same for all five energies. The contibution of this term 
to the total yield is estimated to be about 10%. The nu- 
merical values for dn/dy, T and (nit) — m are listed in 
Table III. 

TABLE III. Particle yield dn/dy, inverse slope T and 
mean transverse mass (nit) — fn for antiprotons and pro- 
tons in central Pb+Pb collisions at 20, 30, 40, 80 and 158 
A GeV beam energies. The errors are statistical. 



Ebeam dn/dy T (m t ) - m 

(AGeV) (MeV) (MeV/c 2 ) 



P 



P 



158 


1.66 + 0.17 


291 ± 15 


384 


± 


19 


80 


0.87+0.07 


283 ± 30 


385 


± 


41 


40 


0.32 + 0.03 


246 ± 35 


355 


± 


51 


30 


0.16 + 0.02 


290 ± 45 


395 


± 


60 


20 


0.06 + 0.01 


279 ± 64 


394 


± 


60 


158 


29.6 + 0.9 


308 + 9 


413 


± 


13 


80 


30.1 + 1.0 


260+ 11 


364 


± 


16 


40 


41.3+1.1 


257+ 11 


367 


± 


16 


30 


42.1 + 2.0 


265+ 10 


362 


± 


14 


20 


46.1 + 2.1 


249 + 9 


352 


± 


13 



DISCUSSION 

The midrapidity mt-spectra for 158 A GeV Pb+Pb col- 
lisions become progressively flatter from peripheral to 
central collisions (Fig. 2 and Table II). Both the mean 
transverse mass (m t ) — m and the inverse slope param- 
eter T increase towards central collisions, although this 
trend is slightly different for p and p, namely the proton 
inverse slope increases somewhat faster, as clearly seen 
in Fig. 4. 

A deviation from a single exponential at low m t values, 
the so called "shoulder- arm shape", is most pronounced 
in central Pb+Pb collisions (Fig. 2) and observed at all 
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FIG. 4: (Color online) The mean transverse mass (rn t ) — rn 
for p and p as a function of the number of wounded nucleons 
N woun d at midrapidity in 158A GeV Pb+Pb collisions. The 
errors are statistical. 



five energies (Fig. 3). These features are generally ex- 
plained by a strong collective ex pan sion (radial flow) in 
hydrodynamical models 0, |U l33l | . A blast wave pa- 
rameterization based on these models provides in fact as 
good a description of the mj-spectra |34j ] as the double 
exponential function. Fitted parameters yield an average 
radial flow velocity of about 0.5 c and a temperature of 
120-140 MeV in the expanding fireball at its disintegra- 
tion stage. 

Interestingly, the p and p transverse mass spectra re- 
veal within errors quite similar shapes in spite of the fact 
that at these energies a significant fraction of protons 
carry baryon number from the incoming nuclei, while 
antiprotons are predominantly pair produced. This ob- 
servation may indicate similar expansion dynamics for 
both particles comprised of thermal motion and collec- 
tive radial flow that determine the transverse spectra of 
particles at the late stage of the expansion. As discussed 
in the introduction, substantial annihilation of primor- 
dially produced p might occur in the baryon rich fireball. 
This effect is expected to lead to a decrease of the p/p 
ratio at uit — m below about 0.4 GeV/c 2 |l5lll7| . In par- 
ticular, at 158A GeV the close similarity of the shapes 
of the measured p and p spectra do not support such a 
prediction. At lower energies a decrease of the order of 
30% cannot be ruled out due to the larger uncertainties 
of the measurements. 

In contrast to the strong centrality dependence, the 
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FIG. 5: (Color online) The mean transverse mass (m t ) — rn 
for p and p at midrapidity as a function of the center-of-mass 



energy per nucleon pair 



' in central Pb+Pb collisions at 



SPS energies (NA49) together with data for central Au+Au 
collisions from the AGS and RHIC HI- The errors 

are statistical. 



shape of the measured transverse mass spectra does not 
change noticeably with beam energy at SPS energies. 
This is illustrated in Fig. 5 where the (mt) — m values 
for p and p in central Pb+Pb collisions at the five SPS 
energies are plotted versus the nucleon-nucleon center-of- 
mass energy. Also shown are measurements at lower and 
higher energies from the AGS [HSI^ and RHIC [H, 
respectively. It is seen that the mean transverse mass 
for both p and p steeply rises with energy at the AGS, 
remains approximately constant at the SPS and resumes 
a slow rise towards RHIC energies. A similar non trivial 
energy dependence was also observed by NA49 for kaon 
and pion production 30] which was considered a possible 
manifestation of the coexistence of partonic and hadronic 
phases or more generall y a very soft equation of state in 
the SPS energy regime |39l l4fJ, |41| . The behavior of the 
m t spectra may thus be taken to support the indication 
for an onset of deconfinement in the early stage of the col- 
lisions at low SPS energies as obtained from the energy 
dependence of pion and kaon yields [Tij . 

The measured particle yields dn/dy differ significantly 
between p and p. Both increase with increasing central- 
ity, but the proton yield rises about two times faster. The 
yields of p and p per wounded nucleon (dn/ dy) / (N woun d) 
as a function of N woun d for the six centrality samples are 
depicted in Fig. 6. An increase of this ratio with cen- 
trality for protons is seen which can be understood as 
a consequence of the increased baryon stopping in cen- 
tral Pb+Pb collisions [H HJ HI]. On the other hand, 
the p yield per wounded nucleon shows a decrease with 
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FIG. 6: (Color online) The yield dn/dy of p and p per number 
of wounded nucleons N woun d and the p/p ratio plotted as 
a function of N woun d at midrapidity in 158A GeV Pb+Pb 
collisions. The errors are statistical. 



increasing centrality. This could be a sign of increasing 
antiproton absorption in central collisions. This effect 
was studied with the microscopic models UrQMD [17] 
and HSD [5] which comprise string excitation, hadroni- 
sation and rescattering (transport) stages. The UrQMD 
results [17] suggest a concurrence of the p enhancement 
in Pb+Pb collisions with respect to p+p reactions and 
an increasing antibaryon absorption with centrality due 
to the growth of the net baryon density at midrapidity. 
However, the observation of nearly identical spectra for 
p and p does not support the UrQMD model calculation. 
The HSD model explains the weak centrality dependence 
of the p yield per wounded nucleon by the contribution of 
multi- meson reactions in the transport phase [5]. These 
regenerate p and thus counteract the loss by absorption. 

The different dependence of p and p yields on N wounc i 
results in significant centrality dependence of the p/p ra- 
tio as shown in the bottom panel of Fig. 6. This ratio 
steadily decreases with centrality by a factor of about 
two within the considered centrality range. A similar re- 
sult was obtained in Au+Au collisions at 11.7^4 GeV |37| . 
The data for smaller collision systems indicate an even 
larger pbar/p ratio than observed in peripheral Pb+Pb 
reactions reaching values of nearly 0.3 for p+p ^2] and 
p+Be, p+S and p+Pb 0] interactions. In contrast 
the p/p ratio measured at RHIC shows very weak vari- 
ation with centrality at the collision energy ^/s NJV = 130 
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FIG. 7: (Color online) The p/p ratio at midrapidity as a 
function of the center-of-mass energy per nucleon pair ^/s NN 
in central Pb+Pb collisions at SPS energies (NA49) together 
with the data for lower e nergies at the AGS p37| and hi gher 
energies at the RHIC |44l l45l |4?| . respectively. The data |4H| 
at top SPS energy are also shown. 



GeV |4J, |45| and almost no dependence on centrality at 
ys^7=200 GeV 0,143. Thus, the net baryon density, 
which is significant at lower energies, strongly affects the 
midrapidity p and p abundances through baryon stopping 
and baryon-antibaryon annihilation in central collisions 
at the AGS and SPS, but not at RHIC energies. 

Fig. 7 displays the midrapidity p/p ratio for central 
collisions as a function of energy at AGS, SPS and RHIC. 
The numerical values are listed in Table IV. The p/p 
ratio rises steeply within the SPS energy range by nearly 
two orders of magnitude. The figure illustrates how the 
collisions evolve from producing a net baryon-rich system 
at the AGS through the SPS energy range to an almost 
net baryon-free midrapidity region at the RHIC. 

Measurements of A production in central Pb+Pb col- 
lisions |U H3 make it possible to analyse strange and 
nonstrange antibaryon production for all five beam en- 
ergies. Of particular interest is the A/p ratio which was 
briefly discussed in the introduction. Note, that the A 
yields used for calculations contain the contribution from 
electromagnetic decays of S° hyperons, which are experi- 
mentally indistinguishable from those created in primary 
interactions. 

The measured values of the A/p ratio for central 
Pb+Pb collisions at 20, 30, 40, 80 and 158 A GeV are 
listed in Table IV and plotted in Fig. 8 together with 
those from AGS and RHIC. The AGS experiments re- 
ported s. A/p ratio of about 3-3.5 for central Au+Au 0,0 
and Si+Au [j| collisions at beam momenta of 11. 7 A and 
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FIG. 8: (Color online) The A/p ratio at midrapidity as a 
function of the center-of-mass energy per nucleon pair ^/s NN 
in central Pb+Pb collisions at SPS energies (NA49) together 
with the data from AGS 0, @ and RHIC ,1.1]. The total 
error bars are drawn. The curve shows the prediction of the 
statistical hadron gas model |4flj . 



14.6^4 GeV/c, respectively. As illustrated in Fig. 8 the 
measurements at the SPS indicate a gradual increase of 
the A/p ratio from 158 A GeV to 30 and 20 A GeV, and 
tend to corroborate the large values for this ratio found 
at AGS energies. At 158^4 GeV the published predic- 
tion 0| for the midrapidity A/p ratio from the UrQMD 
model, which takes into account antibaryon absorption, 
agrees well with the measured value. Predictions for the 
full energy range are not yet available in the literature. 
Since both A and p are newly produced baryons having 
no valence quarks in common with the projectile nucle- 
ons we may compare the midrapidity ratio with the full 
phase space multiplicity ratio predicted by the statistical 
hadron gas model |4flj which uses a smooth parameter- 
isation of the energy dependence of the baryochcmical 
potential. As demonstrated by the curve in Fig. 8 the 
hadron gas model underpredicts the ratio but shows a 
rise towards lower energies similar to the measurements. 
Similar predictions were obtained within non-equilibrium 
versions of the hadron gas model [50L l5l| . 

The increase of the A/p ratio towards top AGS en- 
ergy may also find an explanation in a quark coalescence 
model scenario. Note first of all that we deal here with 
an antihyperon to antiproton maximum because this ra- 
tio has to fall down again at yet lower energies (where 
no data exist due to insufficient statistics) owing to the 
higher A production threshold. Such an antihyperon 
maximum is reminiscent of the maximum in the same 
energy range that was reported recently for the K + /tt + 



FIG. 9: (Color online) Impact parameter dependence of the 
A/p ratio at midrapidity in central Pb+Pb collisions at 158^4 
GeV. The errors are statistical. The curved line pres ents the 
result of hadronic cascade (UrQMD) calculation 18\. 



production ratio 19, 20]. The steep maximum in the rel- 
ative strangeness production was predicted as a signal of 
the onset of deconfinement j^] . In fact if hadronization 
occurs by quark coalescence at the QGP phase boundary 
the A/p ratio essentially reflects the ratio of s to u quark 
densities, while the K + /n + ratio follows from the s to d 
ratio. The d quark density is expected to be proportional, 
in turn, to the u quark density. These considerations thus 
provide a possible explanation of the similar rise of the 
A/p and K + /n + ratios with decreasing collision energy. 

Table IV. The p/p and A/p ratios at midrapidity in cen- 
tral Pb+Pb collisions at SPS energies. The errors are 
statistical. Preliminary results for the A yields at 20 and 
30 A GeV El are used. 



E beam {AGeV) 


p/p ratio 


A/p ratio 


158 


0.058 ±0.005 


1.09 + 0.15 


80 


0.028 + 0.003 


1.22 + 0.14 


40 


0.0078 + 0.0010 


1.31 + 0.19 


30 


0.0038 ±0.0008 


1.81 + 0.37 


20 


0.0013 + 0.0002 


1.72 + 0.58 



For 158A GeV Pb+Pb collisions the data allow a study 
of the centrality dependence of the midrapidity A/p ratio 
which has been suggested to be particulary sensitive to 
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the interplay between production and subsequent absorp- 
tion. As data on the centrality dependence of A produc- 
tion are not yet available from NA49, numerical values 
for the ratio were derived from the yield of A in the most 
central bin 1 measured by NA49 [22j and the calculated 
yields for each centrality bin 2-6 using the NA57 result 
[23j that the A yield within the errors is proportional to 
the number of wounded nucleons N woun d- As shown in 
Fig. 9 the A/p ratio increases from peripheral to central 
collisions by about a factor 2. A much stronger increase 
with centrality has been observed at the AGS 0, • In- 
terestingly, the results at 158^4 GeV are in agreement 
with a UrQMD model calculation which incorporates 
antibaryon absorption. 

SUMMARY 

The midrapidity transverse mass distributions for an- 
tiprotons were measured in Pb+Pb collisions at 20, 30, 
40, 80 and 158 A GeV filling the gap in the data between 
AGS and top SPS energies. They are compared to the 
relevant data for proton and antilambda production. 

The shapes of m t -spectra show a strong dependence 
on collision centrality. In central collisions and at small 
m t the spectra exhibit a pronounced deviation from 
the exponential Boltzmann shape. The data reveal the 
shoulder-arm structure characteristic of prominent radial 
collective expansion flow. The observation of nearly iden- 
tical spectral shape for antiprotons and protons indicates 
a similar expansion history for both particle species. 

No visible change in the mean transverse mass for p 
and p was found in central Pb+Pb collisions within the 
measured SPS energy range. This is similar to what has 
been observed in kaon and pion production and could 
be attributed to the possible formation of a deconfincd 
phase in the early stage of the collision. 

The yield of p normalized to the number of wounded 
nucleons increases with centrality while this quantity for 
p exhibits a decrease. The former effect can be under- 
stood as the result of an increase of baryon stopping with 
collision centrality which leads to contraction of the pro- 
ton rapidity distribution around midrapidity. The de- 
crease of the yield of p per wounded nucleon may point 
to some contribution from annihilation. 

The p/p ratio is found to increase by almost two orders 
of magnitude in central Pb+Pb collisions at 158A GeV as 
compared to 20A GeV beam energy, reflecting the rapid 
decrease of the net baryon density at midrapidity with 
collision energy. 

In central Pb+Pb collisions the A/p ratio shows a 
steady increase with decreasing beam energy approach- 
ing a value of almost 2 at the energies of 20 and 30 A 
GeV. This confirms earlier evidence for this ratio to sig- 
nificantly exceed unity at AGS energies. There may be 
an analogy to the K + maximum also observed in the 



domain of top AGS to lowest SPS energies. 

No fully satisfactory theoretical description of an- 
tibaryon production is available in the literature. The 
new experimental results on p production should help 
to elucidate the possible role of annihilation processes in 
antihyperon production in collisions of heavy nuclei. 
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